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ABSTRACT 

We study and analyse low-resolution spectra of the unordinary late WN star FSZ35 in 
M33. We classify the object as a hydrogen-rich WN8 star. Using the radiative transfer 
code CMFGEN, we determine the physical parameters of this object and compare 
them to the parameters of other WN8 stars including the LBV star Romano's star 
during the minimum of brightness. Unlike Romano's star, the object is fairly stable 
both spectrally and photometrically, that may be attributed to its more advanced 
evolutionary stage or lower luminosity. FSZ35 is shown to possess a compact nebula 
producing faint but detectable [O iii] emission. Location of this object at a large dis- 
tance lOOpc) from the nearest association suggests the object may be one more 
example of a massive runaway star. 

Key words: galaxies: individual: M33 - stars: Wolf-Rayet - stars: supergiants - 
stars: individual: FSZ35 (M33) 



1 INTRODUCTION 

Since middle 1990s, radiation transfer codes designed for 
modeling extended expanding atmospheres became a reli- 
able instrument for studying hot stars with high mass-loss 
rates. These codes make it possible to calculate the at- 
mospheres of outstanding objects like O-supergiants, Wolf- 
Rayet (WR) stars and even Luminous Blue Variables (LBV). 
In particular, the CMFGEN (Hillier & Miller 1998) code 
was used to model the spectra of t] Car (Groh et al. 2010), 
AG Car (Groh et al. 2009) and AFGL2298 (Clark et al. 
2009). CMFGEN and PoWR (Hamann & Grafener 2003) 
were used to determine the properties of the majority of the 
Galactic nitrogen-rich Wolf-Rayet (WN) stars (Crowther et 
al. 1995; Hamann et al. 2006) including WN stars near the 
Galactic center (Martins et al. 2007; Hamann et al. 2010). 

High-luminosity stars are rare objects, and their studies 
in our Galaxy are complicated by dust absorption, crowded 
stellar fields and unknown distances. Therefore, it is impor- 
tant, and sometimes more convenient, to study massive stars 
in nearby galaxies. Studying WR stars that are probable 
progenitors of Type Ibc Supernovae and gamma-ray bursts 
is of especial importance (see review by Crowther (2006) 
and references therein). It is even more difficult to deter- 
mine the distances toward Galactic WN8 stars because they 
can not be associated with clusters or associations (Crowther 
et al. 1995b; Marchenko et al. 1998). One possible scenario 
proposed to explain the apparently enigmatic characteris- 
tics of WN8 stars is that they are runaways, either ejected 
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via dynamical interaction from the cores of forming dense 
star clusters (e.g. Evans et al. (2010)), or accelerated by a 
slingshot-type mechanism during the supernova explosion of 
the initial primary component in a binary system (De Bon- 
der et al. 1997). 

The comprehensive catalogue of Massey & Johnson 
(1998) lists in total 117 nitrogen-rich Wolf-Rayet stars in 
M33 (including WN candidates) . But models of atmospheres 
are constructed for only few objects. Parameters of the 
atmospheres of MCAl-B and B517 were determined with 
the iterative technique of Hillier (1990) (precursor of CM- 
FGEN) in Smith et al. (1995) and Crowther et al. (1997), 
respectively. Using the WM-BASIC code (Pauldrach et al. 
2001), Bianchi et al. (2004) used the ultraviolet spectra of 
six known late WN (WNL) stars in M33 to estimate their 
fundamental parameters: wind terminal velocities, luminosi- 
ties, effective temperatures, radii, mass- loss rates and C/N 
abundance ratios. In the current work, we use the CMF- 
GEN code (Hillier & Miller 1998) to study the spectrum of 
the little-studied WN star FSZ35 in M33. 

The object FSZ35 (a = 0l''33™00".20, 5 = 
-|-30°30' 15". 2, (J2000 epoch)) is situated near the associ- 
ation 128 OB'^ in the M33 galaxy 14' away from galactic 
centre. Projected distance toward the association is 30-40" 
or 100-120 pc. Ivanov et al. were the first to obtain pho- 
tometrical data for this object (IFM-B 174, in their nota- 
tion) in the U, B, and V bands and publish an identifica- 



^ From the catalogue of OB associations by Humphreys 
Sandage (1980) 
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Table 1. Observational log for the SCORPIO data. S/N is signal-to-noise ratio, PA is position angle. 



Exposure Spectral Spectral 
Date time Grism range SX S/N Seeing standard PA 
W [A] [A] ["] star [°] 

FSZ35 

4. 10. 2007 900 -I- 1200 VPHG1200G 4000-5700 5.5 19 1.2 G191B2B -182 

V532 

5. 10. 2007 3 X 900 VPHG1200G 4000-5700 5.5 34 1.1 BD25d4655 -141 



tion chart for it (Ivanov et al. 1993). The star is listed in 
the Ha emission-line object catalogue (Fabrika, Sholukhova 
& Zakharova 1997) as object number 35. In 1998, Massey 
& Johnson published a list of 22 Wolf-Rayet stars selected 
by photometry in three non-standard filters: WN (A4686), 
WC (A4650) and CT (A4752, continuum). Identification was 
then confirmed spectroscopically (Massey & Johnson 1998). 
In this work, the star is listed £is object El and is classified 
as a Wolf-Rayet star of nitrogen sequence, WN8 subtype. 
Sholukhova et al. (1999) note similarity between the spectra 
of FSZ35 and V532 (Romano's star). The latter is a better 
studied object located in the outer spiral arm of M33 at a 
distance of about 17' from the centre. V532 is universally 
recognized as an LBV star. It demonstrates pronounced pho- 
tometrical and spectral variability (Kurtev et al. 2001; Pol- 
caro et al. 2003; Maryeva & Abolmasov 2010). Hence we pay 
especial attention to comparison between these two objects. 
Similar modeling for V532 will be published in a separate 
paper (Maryeva & Abolmasov 2011). 

This paper is organized as follows. The observational 
data and data reduction process arc described in the next 
section. We characterise and classify the observed spectrum 
in Section 3. In Section 4 we describe the basic properties 
of the CMFGEN code and its main assumptions, while in 
Section 5 we present and analyse the modeling results. In 
Section 6 we discuss the results. Finally, in Section 7 we 
summarize the main points of our work. 



2 OBSERVATIONS AND DATA REDUCTION 

In this work, we use a spectrum of FSZ35 obtained at the 
6m Special Astrophysical Observatory (SAO) telescope^. 
Two exposures, 1200s and 900s in length, were obtained 
with the SCORPIO multi-mode focal reducer (Afanasicv & 
Moiseev 2005) in the long-slit mode on October 4, 2007. 
VPHG1200 G grism was used providing the spectral range 
of 3950-5500 A. The data were reduced using the ScoRe 
package for long-slit data reduction, written in IDL espe- 
cially for SCORPIO long-slit data reduction. This package 
consists of procedures created by V.Afanasiev, A. Moiseev, 
P.Abolmasov and O. Maryeva. Package includes all the stan- 
dard stages of long-slit data reduction process. The final 

^ Spectral data were taken from the archive of Special Astro- 
physical Observatory (SAO) of Russian Academy of Sciences, 
http://www.sao.ru/oasis 



spectrum has spectral resolution ~ 5.5 A (weakly depen- 
dent on wavelength) and signal-to-noise ratio per resolution 
element in continuum ~ 20. 

Besides the spectrum of FSZ35, in this work we use a 
spectrum of Romano's star obtained at the 6m SAO tele- 
scope practically at the same time with the same grating 
under similar conditions. Details about the spectral data on 
V532 and data reduction may be found in our recent paper 
Maryeva & Abolmasov (2010). Basic information about the 
observational data on both objects is summarized in Table 1. 



3 SPECTRAL CLASSIFICATION 

Figure 1 shows the spectra of V532 and FSZ35 obtained 
under similar conditions. The spectral appearance of FSZ35 
shows strong similarities with V532. The spectrum of FSZ35 
is as rich as the spectrum of V532 in emission lines, but 
the presence of nebular lines is questionable (see below Sec- 
tion 6) and hydrogen lines are fainter. Table 2 presents the 
lines detected in the spectrum of FSZ35. Lines of helium 
and hydrogen in the spectrum of FSZ35 are much broader 
in comparison with the these lines in the spectrum of V532. 
For example, FWHM (Full Width at Half Maximum) of H/3 
is 12.0±0.2 and 5.5±0.lA for FSZ35 and V532, respectively. 
For HelA4921, FWHM arc 9.5 ± 0.2 and 4.4 ± O.lA. 

Using nitrogen and He II lines, we classified V532 in Oc- 
tober 2007 as a WN8 star using the classification of Smith, 
Crowther & Prinja (1994) for WN6-11 stars. Locations of 
V532 on the diagrams "equivalent width of He I A5876 ver- 
sus He II A4686" and "equivalent width of He II A4686 versus 
FWHM of this fine" are consistent with the WN8 subtype 
(Maryeva & Abolmasov 2010; Polcaro et al. 2010). 

Similar consideration allows to classify FSZ35 as a WN8 
as well. Besides this, the NlvA4057 line clearly seen in our 
spectrum is never present in WN9 spectra, thus exclud- 
ing FSZ35 identification as a later-subclass object. Unfor- 
tunately, our data cover only a limited spectral range from 
3900 to 5500A. Therefore we can not use the Hex A5876 fine 
for spectral classification. 

We used a quantitative chemistry-independent criterion 
based on the FWHM of the He ii A4686 Une for alternative 
spectral classification (see, for example. Smith et al. (1995)). 
Equivalent width (EW) of the He II A4686 line in the spec- 
trum of FSZ35 is 20A. FWHM of this line is 10.5 ± 0.5A. 
In Fig. 2, we show the location of V532 and FSZ35 on the 
diagram of the EW of Hell A4686 versus the FWHM of this 
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Figure 1. Spectra of FSZ35 (solid line) and V532 (dotted line) obtained in October 2007. Tlie spectra are normalized by the local 
continuum level. 

Table 2. List of emission lines detected in the spectrum of FSZ35. For higher signal-to-noise ratios, equivalent widths (EW) are given. 



A, A 


Ion 


EW, A 


A, A 


Ion 


EW, A 


3964.73 


Hex 




4523.60 


NIII 




3970.08 


He-)-HeII 




4530.80 


NIII 




3994.99 


Nil 




4534.60 


NIII 




4009.00 


He I 




4541.60 


Hell 




4025.60 


He I-I-He II 




4547.30 


NIII 




4057.80 


NIV 


1.10 ±0.15 


4601.50 


Nil 




4088.90 


Si IV 


1.2 ± 0.3 


4607.20 


Nil 




4097.31 


NIII 




4613.90 


Nil 




4101.74 


H(5+HeII 


1 4.4 ±0.4 


4621.40 Nil 




4103.40 


NIII 




4630.54 Nil 




4116.10 


Si IV 


2.4 ±0.4 


4634.00 


NIII 


7.2 ±0.1 


4120.99 


He I 




4640.64 


NIII 


4.9 ±0.1 


4143.76 


He I 


1.0 ±0.3 


4643.09 


Nil 




4199.80 


Hell 




4650.16 


C HI 


3.3 ±0.1 


4236.93 


Nil 




4685.81 


He II 


20.0 ±0.5 


4241.79 


Nil 




4713.26 


He I 


2.8 ±0.2 


4241.79 


Nil 




4861.33 


H/3 


13.4 ±0.3 


4340.47 


H7-|-He II 


5.5 ± 0.4 


4921.93 


He I 


4.1 ±0.2 


4387.93 


He I 




5001H-5007 NlI+[OlIl] 




4471.69 


He I 


3.0 ±0.3 


5015.67 


He I 


2.9 ±0.2 


4510.90 " 






5314 " 






4514.90 


NIII 


~ 1 


5320 


NIII 




4518.20 






5327 
5411.50 


Hell 


1.8 ±0.2 
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Figure 2. The location of FSZ35 and V532 (October 2007) on 
equivalent width versus FWHM diagram for the He nA4686 line. 
Known Galactic (filled symbols) and LMC (open) WN stars are 
shown for comparison: WN7 by squares, WN8 by circles, WN9 
by diamonds, and WNIO by triangles. Data on these objects were 
taken from Crowther & Smith (1997); Crowther & Bohannan 
(1997). 



line. Position of the object in this diagram is consistent with 
our identification of FSZ35 as a WN8 star. 

We confirm the spectral classification of the object as a 
WN8 star. This our result also suggests the object is rrmcli 
less variable than V532 that spans a broad range of effective 
temperatures and luminosities and becomes a late WN star 
only in deep minima. FSZ35 is thus closer to "quiet" WN 
stars. 



4 MODELING 

For our analysis we used the non-LTE radiative transfer code 
CMFGEN (HilHer & Miller 1998). CMFGEN solves radia- 
tive transfer equation for objects with spherically-symmetric 
extended outflows using either the Sobolev approximation 
or the full comoving-frame solution of the radiative transfer 
equation. To facilitate simultaneous solution of the transfer 
equation and statistical equilibrium equations, partial lin- 
earization method is used. To facilitate the inclusion of metal 
line blanketing in CMFGEN, superlevel approach (Anderson 
1989, 1991) is used. In this formalism, levels with similar 
properties are treated as one and have identical departure 
coefficients. This allows to save considerable amount of com- 
puter memory and time. Recent versions of the code incor- 
porate also the effect of level dissolution, influence of reso- 
nances on the photoionization cross section, and the effect 
of Auger ionization. 

Clumping is incorporated into CMFGEN using volume 
flUing factor approach (Hillier & Miller 1999). Filling factor 
is allowed to depend on radius. By default, the wind is con- 
sidered homogeneous at the hydrostatic radius and becomes 
more and more clumped with the wind velocity. Clumping 
tends to reduce the derived mass-loss rates by a factor of 
~ 3 — 5 (Marcolino et al. 2007). Unclumped mass loss rate 
(i. e. calculated not taking clumping into account) is related 



to the volume filling factor / by relation Mund = Mc\/^/f 
(see for example Herald et al. (2001)). 

Every model is defined by a hydrostatic stellar radius 
R,, luminosity L*, mass-loss rate Md, filling factor /, wind 
terminal velocity Voo , stellar mass M, and by the abundances 
Zi of included elementary species. We assumed a constant 
turbulent velocity of 20 km/s. Its variations do not affect 
the resulting spectrum mtich save for the equivalent width 
of HeIA4686 that becomes slightly brighter in more turbulent 
atmospheres. 



5 MODELING RESULTS 

Using the model of Romano's star calculated by Maryeva & 
Abolmasov (2011) as the seed model we adjusted its param- 
eters to reproduce the observed spectrum of FSZ35. As for 
V532, H, He, C, N, Si, Fe, O, Ne, Mg, S, Ar, Ca, and Na were 
included in calculations. We considered only sub-solar iron 
abundances (between 0.2 and 0.5 solar). Hydrogen lines in 
the spectrum of FSZ35 are weaker than in the spectrum of 
V532, suggesting lower H/He values of ~ 0.6^ 0.8 (by num- 
ber), unlike V532 where hydrogen abundance is estimated 
as H/He~1.3. 

We increased the effective temperature of the model in 

order to reproduce the NlvA4057 emission present in the 
spectrum, and then varied the luminosity at a constant ef- 
fective temperature level to fit the observed V-band lumi- 
nosity. V = 18"". 7 is reported by (Massey & Johnson 1998), 
in consistence with the observed spectrum slightly hotter 
than that of V532 in deep minimum. In our calculations we 
supposed that distance toward M33 is _D = 847 ±60 kpc and 
the distance modulus is (m - M) = 24.64 ± 0.15™ (Galleti 
et al. 2004). We neglected intrinsic extinction of M33, since 
the object is distant from center of the galaxy. Extinction 
in the Galaxy is estimated as E{B — V) = 0™. 029 (corre- 
sponding to Ay = 0?1 for Rv = 3.1) by the NED extinction 
calculator (Schlegel et al. 1998). 

At the next step, we change the mass loss rate and 
terminal velocity Voo- Lower resolution and spectral range 
do not allow to make reliable estimates of the terminal wind 
velocity, as we did for V532 in Maryeva & Abolmasov (2010). 
Terminal velocity given in Table 3 was found by fitting the 
observed line profiles. The filling factor was set to / = 0.1 
for all the models. At the last step of the modeling process, 
we adjust the nitrogen abundance. 

In Figure 3 we present the spectrum of FSZ35 and 
our model. We succeed in achieving good agreement with 
the observational equivalent widths and profiles of hydrogen 
linos and singlet lines of neutral helium as well as the line 
He IA4713. 

In the best-fit model, H/He ratio is 0.8 and, conse- 
quently, mass fraction of hydrogen is ~0.17. FSZ35 con- 
forms to the definition of H-rich WN stars (Nugis & Lamors 
2000), that must have hydrogen mass fraction > 0.05. The 
mass loss rate is Md = 2.6 • lO~^M0/year and Mund = 
8.22 • 10~^Mo/year. The mass loss rate is given in Table 3 
together with other wind and stellar parameters. For com- 
parison, the values of these parameters for some other WN 
stars taken from the literature are given in the table. Note 
that the parameters of WR116 and WR89 were calculated 
using PoWR code (Hamann & Grafener 2003) while the pa- 
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Figure 3. The spectrum of FSZ35 (solid line) compared with the best-fit CMFGEN model (dotted line). 



rameters of the other stars using CMFGEN. Table 3 shows 
that the parameters of FSZ35 arc similar to these for some 
other WN8 stars given in the table, irrespective of a code 
which was used the estimate the parameters. In our calcula- 
tions wc neglected intrinsic extinction in M33. If extinction 
near FSZ35 is A,, = 0™. 1, luminosity increases by about 10% 
and becomes L, — 6.3 • IO^Lq. 

FSZ35 differs from V532 in the values of mass loss 
rate (mass-loss rate of V532 is M = 1.6 • lO~®M0/year). 
This difference is quite normal, because mass-loss rate in- 
creases with luminosity and with helium abundance (Nugis 
& Lamers 2000). Moreover the nitrogen abundance N/He is 
IQ-^ (~ 2.7 solar) unlike V532, where N/He is 3.0 x 10"^ 
(~ 7.1 solar). 



6 DISCUSSION 

6.1 Nebular Contribution 

Although the spectra of V532 and FSZ35 are very 
similar, a close look at the subtle differences can 
be illuminating. All the forbidden lines, including 
[FeIIl]AA4658.1, 4701.5, 5270.4, arc cither absent or much 
fainter for FSZ35. First we suggested the line near the 
HeiA5015 emission belongs to the NiiAA4987 - 5007 blend 
but we did not find reasonable solution that could ac- 
count for the whole nitrogen spectrum (NllAA5001 — 10, 
N IIIAA4634, 40 and NivA4057). The Une position is consis- 
tent with the [O iii]A5007 emission. Total flux in this line is 



F([OlIl]) ~ (2.8±0.2) X 10"^^crgcm"^s"^ The other com- 
ponent of the doublet, [OlIl]A4959, should thus have the 
flux of ~ 9.3 X 10~^^erg cm~^s~^ and is unreachable in the 
spectrum. 

Using the best-fit CMFGEN model spectrum as input, 
we run a set of Cloudy (Ferland et al. 1988) version 08.00 
photoionization models. We restrict the size of the puta- 
tive nebula around FSZ35 with the outer radius of 0.5 pc 
(there are no indications for resolved diffuse nebular emis- 
sion around the star) and fix the inner radius to 0.1 pc. 
Metallicity pattern Hii region was used with the metal- 
licity decreased by a factor of 2. We do not change nitrogen 
and helium abundances in the nebula because they can be 
less affected by chemical evolution of the star and represent 
rather the chemical composition of its early hydrogen-rich 
ejecta. 

The only varied parameter was hydrogen density. For 
several values of this parameter, we calculate the luminosi- 
ties of several nebular lines and give them in figure 4. For 
comparison, the observational estimate for the luminosity 
in [Oiii]A5007 is given. Evidently, the observed value is in 
good agreement with the predicted nebular luminosity for 
hydrogen density tih ^ SOOcm^''. Note however that oxygen 
abundance may be lower in the ejecta consisting the nebula 
(see below). If oxygen is under-abundant by a factor of sev- 
eral, the observed flux is consistent with uh ~ 10■^cm^'^ 

We also expect emission lines of [Sii] and [Nil] in the 
red spectral range to be much fainter. However, they may be 
stronger if the matter in the nebula is denser and chemically 
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Table 3. Derived properties of FSZ35 and V532 in the minimum of brightness (October 2007), and a comparison with similar stars in 
Milky Way (MW) . is mass fra<;tion of hydrogen 



Star 


Gal. 


Sp. 




R* 


logL* 


log Muncl 




Xh 


Ref 






type 


[kK] 


[Rq] 


[Lq] 


[Me/year] 


[km / s] 


[%] 




WR124 


MW 


WN8h 


32.7 


18.0 


5.53 


-4.2 


710 


13 


[1] 


WR116 


MW 


WN8h 


39.8 


21.0 


6.0 


-3.7 


800 


10 


[2] 


WR89 


MW 


WNSh 


39.8 


26.5 


6.2 


-4.2 


1600 


20 


[2] 


WR40 


MW 


WN8h 


45.0 


10.6 


5.61 


-4. 


840 


15 


[3] 


WR16 


MW 


WNSh 


41.7 


12.3 


5.68 


-4.3 


650 


23 


[3] 


FSZ35 


M33 


WNSh 


36.5 


19 


5.76 


-4.1 


SOD 


16.5 




V532 


M33 


WNSh 


34.0 


20.8 


5.7 


-4.3 


360 


24 


[4] 



[1]- Crowther et al. (1999), [2]- Hamann & Grafener (2003), 
[3]- Herald et al. (2001), [4]- Maryeva & Abolmasov (2011) 



101 F 



100 r 




nH, cm ^ 

Figure 4. Dependence of nebuleu: line luminosities on the hy- 
drogen density of the putative nebula around FSZ35. Solid is 
H/3, dotted is [OIII]A5007, dashed is [Sll]A6717 and [Nll]A65S3 
is shown by a dash-dotted line. Horizontal lines show the mea- 
sured luminosity of the [Olll] line (1<t confidence interval). 

evolved, or shock excitation also contributes to the emission 
of the nebula. 

Relying on the similarity of the central stars, one may 
propose that the nebula is similar to Ml-67 observed around 
WR124 (Crowther et al. 1999). Strongest nebular Unes in 
the spectrum of Ml-67 are lower-excitation lines such as 
[Sii]AA6717,3lA, [Nii]AA6548,84A and [Oii]A3727A. 

Another effect making the putative Wolf-Rayet nebula 
of FSZ35 hard to detect with the data on hand may be oxy- 
gen under-abundance (by a factor of 5^10) reported for Ml- 
67 and other Wolf-Rayet nebulae and related objects (see 
Esteban et al. (1991) and references therein). Oxygen deficit 
is a predicted by-product of the CNO cycle and should have 
(as well as carbon and neon deficits) an amplitude similar 
to that of nitrogen over-abundance. 

6.2 Lack of Vciriability 

The spectrum analysed in this work is identical within one 
spectral sub-class to the spectrum published by Massey & 



Johnson (1998). We also classify FSZ35 as WN8. Unlike 
V532, FSZ35 does not demonstrate any prominent spectral 
variability. 

It is also fairly stable photometrically. V-band magni- 
tude of the object was 18™. 7 in 1993 (Massey & Johnson 
1998) and 18™ 8 in 2004 (Hartmarm et al. 2006), that sig- 
nificantly restricts the possible amplitude of photometrical 
variability. In the V band. Am < 0™. 1. The lack of photo- 
metrical variability may be ascribed to the more advanced 
evolutionary status of FSZ35 (Meynet et al. 2011). 

6.3 FSZ35 as a massive runaway 

FSZ35 is located at the distance of about 35" (~ 115 pc) 
from the association OB 128. Suppose that once FSZ35 was 
a member of the association and was ejected via slingshot- 
type dynamical interaction. If its peculiar velocity is ~ 
100km s"\ as for WR124 (Siriarmi et al. 1998; Marcheuko 
et al. 2010), it could have been expelled from the parent 
cluster about a million years ago. 

Offset positions with respect to the probable parent as- 
sociations (at distances ~ lOOpc) and unexpectedly large pe- 
culiar velocities (of the order ~ lOOkm s~^) seem to be com- 
mon for very luminous and massive stars like V532, FSZ35 
and Galactic late WN stars like WR20a and WR124. A sce- 
nario was proposed by Gvaramazde & Gualandris (2011) 
that applies three-body dynamical interaction in the cores 
of young massive clusters and star-forming regions (simi- 
lar to 30 Doradus) to reproduce the observed population of 
massive runaways. This scenario has the disadvantage of re- 
lying heavily on the formation of massive stars in the cores 
of massive young clusters that are, in contrast with LMC, 
absent in M33, where massive stars are rather formed in 
dispersed stellar associations. 

Instead we would rather propose that very massive stars 
are formed in dense groups containing several stars each. 
This is confirmed, for example, by the multiplicity increas- 
ing with stellar mass (Zinnecker & Yorke 2007) both in 
young star clusters and associations. It is reasonable that 
higher fraction of massive binaries will be accompanied by a 
higher fraction of massive rrmltiple systems. When formed, 
such systems are often unstable (Kiseleva et al. 1998) and 
dynamic interaction between its components should both 
produce a larger fraction of runaways at these masses (~ 
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IOOM0) and a larger fraction of binaries. The characteristic 
pecuhar velocities (~ 100km s~^) of these "childhood run- 
aways" may be reproduced if the initial spatial sizes of the 
systems are < lO^^cm. 



7 CONCLUSIONS 

We analyse the optical spectrum of the little-studied WNL 
star FSZ35 in M33. About 40 spectral hncs in the 4000 
55OOA wavelength range are identified. We classify FSZ35 
as a WN8 star, confirming the result of Massey & Johnson 
(1998) and put upper limits for the spectral and photomet- 
rical variability of this object. 

Using non-LTE code CMFGEN wc estimate the param- 
eters of FSZ35 (bolometric luminosity, stellar radius, mass 
loss rate, wind velocity, elementary abundances) and com- 
pare them to the corresponding parameters of other WN8 
stars including the LBV star V532 during the mininmm of 
brightness. FSZ35 is a H-rich WN8 star where the mass 
fraction of hydrogen is 17% (H/He=0.8). The best-fit pa- 
rameters of the model are: luminosity L = 5.3 • IO^Lqi 
mass loss rate 2.1- lO~'M0/year, nitrogen abundance N/Hc 
= 13.5(N/He)0, effective temperature at hydrostatic radius 
T, = 36 470K (7?. = ISEq) and at the Rosseland photo- 
sphere Tij,u=2/3 = 35160K. Derived parameters of FSZ35 
atmosphere correspond to a typical WN8 star. 

We find that FSZ35 has a surrounding nebula, possibly 
of low excitation and deficient in oxygen, that decreases its 
detectability. 

Position of FSZ35 at the outskirts of association 128 OB 
suggests that it was expelled from this association about a 
million years ago at a velocity of ~ 100km s"'^. 
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